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Abstract 

We have revisited the still unresolved puzzle of the dispersion of the Raman 
disordered- induced D band as a function of laser excitation photon energy 
in graphite-like materials. We propose that the D-mode is a combination of 
an optic phonon at the K-point in the Brillioun zone and an acoustic phonon 
whose momentum is determined uniquely by the double resonance condition. 
The fit of the experimental data with the double-resonance model yields the 
reduced effective mass of 0.025me for the electron-hole pairs corresponding to 
the A2 transition, in agreement with other experiments. The model can also 
explain the difference between uss and ujas for D and D* modes, and predicts 
its dependence on the Raman excitation frequency. 
PACS numbers: 63.20.Dj, 63.50.+X, 78.30.Ly 

I. INTRODUCTION 

Carbon based materials, ranging from highly oriented pyrolitic graphite, diamond like 
carbon films, fullerenes and carbon nanotubes have been the subjects of extensive studies 
over the past two decades because of their numerous technological applications. Raman 
spectroscopy is a powerful non-destructive probe for characterising different forms of carbon 
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materials-diamond, crystalline graphite, glassy carbon, hydrogenated amorphous carbon, 
fuUerenes and nanotubes. Raman spectra of these materials gives considerable information 
on the nature of chemical bonding (sp^ vis-a-vis sp^) and relative abundance of sp^ and sp^ 
bonds, residual stresses, crystallite size and dopants. There is a long-standing puzzle in the 
Raman spectra of graphitic materials. Single crystals of pristine graphite {D^f^ space group 
symmetry) has two in-plane Raman active modes: one at 42 cm~^ and the other at ~ 
1582 cm~^ (G-band)0. In presence of disorder, an additional Raman band (called D-band) 
is observed at ~ 1350 cm~^ for a laser excitation energy of 2.41 eV (wavelength 514.5 nm), 
irrespective of the kind of carbon materialSHl. 

The important features of the Raman spectra associated with the disorder- induced D- 
band are as follows: (i) The E^-dependence of the D-band frequency is essentially inde- 
pendent of the type of carbon involved. It occurs in all sp^- hybridized disordered carbon 
materials, powdered and randomly oriented crystalline graphite, glassy carbon, boron-doped 
highly ordered pyrolytic graphite, carbon black, multicomponent carbon films, hydrogenated 
amorphous carbon and more recently, in carbon nanotubesi. (ii) The frequency of the D- 
band shifts upward with increasing excitation laser energjil ■ ujd = 1279 cm~^ for = 
1.17 eV and ud ~ 1410 cm~^ for E^ = 3.54 eV (~ 40 cm~^/eV). Fig.l shows the data from 
various measurements made on various kinds of disordered carboniS and carbon nanotube^ . 

(iii) The frequency of the associated second order (D*) band (~ 2700 cm~^ for Ei=2.41 eV) 
also shifts up with increasing E^,. The dispersion duj£,*/dEL is almost twice that of dun/dEL. 

(iv) The E^-dependant second order D* band is observed in Raman spectra of single crystal 
graphite even though the first order disorder induced D-mode is absenti. (v) The D-band 
intensity decreases as E^ increases from 1.16 eV to 4 eV. In contrast, the intensity of the 
G-band is maximum at E^ ~5 evi. The D band is not observed with excitation energies 
El > 4.2 eV as shown in Fig.2 which displays Raman spectra of graphite powder (with a 
small amount of sodium sulphate for internal intensity calibration) using El = 2.54 eV and 
4.82 eV. The lines marked by asteriks are from internal vibrational modes of sulphate ions. 
The intensities have been normalised with respect to the sulphate mode at 980 cm~^. The 



data in the inset taken from Wang et a/.i show the intensity ratio of D band with respect to 
that of G band for 1.16 < < 4 eV. The decrease of this intensity ratio as Ej^ increases has 
been seen by others as well^. In other words, the resonance excitation profile of the D-band 
will show a maximum at El less than 1.16 eV, which is very different from the G-band 
behaviour. These experimental observations are crucial to arrive at the correct explanation 
for the origin of the D-band. (vi) The peak positions of the D and D* bands observed in 
Stokes (us) and anti-Stokes are not the same@: Aus = ^^as — =7 cm~^ for the D-mode 
and ujas — ^s =28 cm^^ for the D* mode when El =2.41 eV. In order to make sure that this 
anomalous difference between Stokes and anti-Stokes peak positions of the D and bands 
is genuine, it was seen in the same experiment that Au ~ for the G-band, as is usually 
expected. 

It is surprising that even after twenty years of experimental observationS, the origin 
and dispersion of the D-band is still not quantitatively understood. Tuinstra and KoenigB 
attributed the D-mode to the Ai^ type mode of a finite graphite crystallite with symmetry 
lower than Dg/j of graphite, which corresponds to a longitudinal acoustic mode for the 
infinite lattice. Nemanich and Solini and Lespande et a/.i have invoked the breakdown of 
wave- vector selection rules due to disorder or finite size of crystals and correlated the D-band 
features with the density of states of graphite. These models, however, cannot explain the 
strong dispersive nature of the D-mode as a function of El- The observation of second-order 
D* band for single crystal graphite which shows no disorder-induced D mode also rules out 
the explanation of the dispersion to crystallite size selective resonance Raman scattering 
process which occurs in polyacetylene. The recent studies0'i explained the dispersion of the 
D-mode following Boronov et a/.0 who attributed the dispersion of the D-band in terms of 
the coupling between electrons and phonons with the same wave-vector near the K-point of 
the Brillouin zone. This model has a serious drawback that it does not justify the reasons 
for the coupling and cannot address the resonance of the D-mode at El < leV, as well as 
the difference in ujas and us for the D and D* modes. 

Almost three years backi, we have proposed disorder-induced double 
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gap of ~ 1 eV in the graphite band structure to result in the dependence of the phonon 
wavevector q and hence the phonon frequency uj on the energy of the exciting laser energy. 
Very recently, Thomsen and Reich^l (TR) have wrongly questioned the existence of such a 
gap of ~ 1 eV in the electronic structure of graphite. They have also used double resonance 
using the electronic linear bands at K-point in the Brillouin zone. The serious difficulty 
with their calculations is that the calculated intensities of the D-mode for different incoming 
photon energies do not decrease as increases. On the contrary, the calculated intensity 
for El = 4 eV is higher than for E^ = 2 eV (see Fig. 3 of Thomsen and Reichlii). This is 
completely opposite to the experiments as shown in Fig. 2. There is another difficulty we 
have with the TR's calculation. As given by their Eq. (4), the magnitude of q increases with 
E/^. For optic branch along FK direction, the phonon frequency is a decreasing function 
of q and hence uj should decrease with increasing E^ which is contrary to the experiments. 

In this paper, we present a model based on double resonance which can address all the 
novel features associated with the D mode. Keeping in mind that the acoustic phonon 
branches are strongly coupled to the high-frequency optic branches at the K-point in 
graphitelii, we propose that the D-mode is a combination of an optic phonon at the K-point 
in the Brillouin zone and an acoustic phonon whose momentum is determined uniquely by 
"double resonance" condition^. In usual second order scattering, if the optic phonon is 
described by a wave- vector qi and the acoustic phonon by a wave- vector g*2; conservation of 
quasi-momentum requires qi + q2 = 0, where g! and q2 can span the entire Brillouin zone. 
Since D-mode is disorder-induced, we suggest that qi is fixed at the K-point of the Brillouin 
zone (maximum in the density of phonon statesil) and the magnitude of q2 is determined 
by the double resonance conditional. The quasi-momentum conservation is achieved by 
disorder-scattering. In a double resonance process, real (not virtual) transition takes place 
between two electronic states with emission of a phonon. Since the electronic states have 
dispersion, phonons of appropriate wavevectors are needed to keep the electronic transition 
real. Double resonance have been seen for optic phonon overtones in GaAjli, 2Fi2 optic 
phonon jli in CU2O and for two acoustic phonons in Ge0. 



II. MODEL 

Following Martin and Falicov0, we will discuss the two-plionon Raman scattering arising 
from an iterated one-phonon scattering in second order0. Near resonance condition, the 
incident photon is absorbed to create an electron-hole pair, with electron occupying the real 
conduction band states and hole in the real valence state. Assuming parabolic dependence 
of the electron-hole pairs in directions perpendicular to KH, the wave- vector of the electron 
or hole k created by a photon of energy El can be written as -E^ = A + where /i is the 
reduced mass ( - = ^ + ^ I , rnt and ml. are the electron and hole effective masses and A 
is the band gap at K-point in the Brillioun zone. We will now address the nature of critical 
point involved in the resonance Raman scattering of D-band. 

We recall that the D-mode intensity increases as Ej;^ is decreased from 4.2 eV upto 1.16 
eV, suggesting that the D-band excitation profile has a maximum below 1.2 eV. We suggest 
that the appropriate band gap associated with the resonance of the D-mode is at ~ 0.8 eV 
where a maximum is seen in the optical reflectivity measurements on graphiteS'El. This 
gap is associated with the energy difference between Ei and E3 bands at the K point of 
the Brillouin zone as shown in the band structure of graphite (Fig. 3a) taken from Ref. [19- 
21]. The wavevector is in the basal plane. The bands shown here are the A-n bands in 
three-dimensional graphite. Recall that a single sheet of graphite has 2-tt bands with linear 
dispersion and Fermi level passes through the point of interaction (K-point) as depicted in 
Fig. 3b. Band structure calculations treating the interlayer coupling as a small perturbation 
find a small gap at K point in the Brillouin zone. The optical transition between Ei and E3 is 
termed as A2 transition. Within the Slonczewski-Weiss-McClure band model@ the energies 
of these optical transitions are mainly determined by the overlap interaction parameter 71 ~ 
0.4 eV for nearest-neighbour atoms on successive layersEll and energy of A2 transition is ~ 
271. Thomsen and ReichEl have therefore, wrongly questioned the existence of a gap ~ 0.8 
eV in the 3D band structure of graphite as suggested Sood et a/.i 

The mechanism we suggest is that the electron in the conduction band (or the hole 
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in the valence band) is scattered by an optic phonon, changing the electron state from k 
to k'. Another accoustic phonon scatters this electron at k' to k", followed by an impurity 
(disorder) scattering from k" to k'" and the electron-hole recombine to produce the scattered 
photon with a frequency shift given hj u = Uop + oJac- As mentioned before, we proposed 
that the optic phonon involved correspond to the maximum in the density of phonon states 
at the K-point in the Brillouin zone (shown by a solid dot in the phonon dispersion of 
graphite in Fig. 3) and the wavevector q of the acoustic phonon is so chosen to satisfy the 
double resonance condition. This Raman process involving disorder mediated two-phonon 
scattering is fifth order in perturbation theory and will have four energy denominators as 
explicitly written by Kauschke et a/.@, out of which two energy denominators will be the 
same. The resultant expression for the Raman matrix element is to be integrated over k, 
k', k" and k'" . The Raman intensity will be maximum when all the three denominators 
resonate simultaneously. It can be shownEUll that the magnitude of q satisfy the condition 
for Stokes Raman scattering, for A + huop < < A + TioJop + huiac-. 
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{EL-A)-^ + {EL-A-hLOop)-^ (1) 



For El > A + hojop + huac, another resonance becomes possible for 
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{El- A- huop) ^ + {El - A - hujop - hujaM))^ (2) 



For a given E^, the two-phonon Raman peak position for Stokes scattering will be given by 



LOop{K) + UJac{qs) (3) 



with qs given by Eq. (1) or (2). For A ~ 0.8 eV and E^, > 1.16 eV, Eq.(2) is the appropriate 
relation to be used. In graphite, there is a peak in density of states around 1230 cm~^ 
corresponding to the optic phonon frequency at the K point in the Brillioun zonJil. We 
therefore, take Uop{K) to be 1230 cm^^ 

Regarding acoustic phonon, it has been seen that the acoustic phonon branch corre- 
sponding to the branch marked T in the phonon dispersion relation of graphite shown in 
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Fig. 4, is strongly coupled to the optical phononS. We take uJadq) = vtaQ where vta is the 
transverse accoustic phonon velocity (= 1.23x10^ m/s). Since cOac uJop, one can simplify 
Eq.(2) and write 



i^aciqs) = QsVtA = '^VtA {^L - ^ ' f^^^op) ^ (4) 

The solid line in Fig.l is a fit to Eq.(3) and (4), with A and fi as adjustable parameters 
and Uop = 1230 cm~^ The values obtained are A = 0.85 ± 0.05 eV and /i = 0.025 me. 
Taking the electron and hole effective mass to be same, m* = 0.05me, which is in excellent 
agreement with the values obtained from other experimental measurements like cyclotron 
resonanceii. We also tried to fit the data in Fig.l using Eq.(4), with Uop also as an adjustable 
parameter along with A and /i. In this case also the fit is very good with parameters, Uop = 
1211 cm^^, A = 0.75 eV and fi = 0.026me. The values of the parameters in both cases being 
in resonable agreement with those obtained from other experiments gives us confidence in 
the double resonance model to understand the origin of D-band. 

Recent Raman measurements on ion-implanted highly oriented pyrolytic graphite by Tan 
et a/.i show that the peak positions in the Stokes and anti-Stokes spectra for the D band 
are different, with uas > ^s- In our model, using the disorder- induced double resonance, 
wavevector of the phonon in anti-Stokes Raman scattering for the D mode will be given by. 
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{El-a + hujop)^ + {El-a + fiujop + nujaM))^ (5) 



in place of Eq.(2). Therefore 
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Ad = UAs -uJs = 2vTA\h^{{EL - A + o^op) ^ - (^l - A - Uop) ^ } (6) 

V a 

For the mode, two optic phonons and two acoustic phonons will be involved. This 
process will not require disorder- induced scattering to conserve momentum. In this case 
Ad* ~ 2A/). Putting in the values for A = 0.85 eV and fi = 0.025me as obtained from 
fit of Fig.l, ujas ~ for D mode = 13 cm^^ (E^ = 2.41 eV) while the difference observed 
experiment allyS is ~ 7 cm~^. In the case of D* mode we find the difference to be 26 cm~^ 



which is very close to the experimentally observed value of 25 cm~^ of Tan et alB. Recent 
Stokes and anti-Stokes measurements by Zhang et a/.il on carbon nanotubes also shows a 
difference in the peak position of the D band in Stokes and anti-Stokes spectra and this 
difference is excitation laser energy dependent, very similar to the predictions of our model. 

In conclusion, we have addressed the long standing problem of the dispersion of the 
Raman D band as a function of laser excitation energy using a simple model based on disorder 
induced two phonon scattering and the double-resonance. The fit of the experimental data 
in this model yields the reduced effective mass of the electron-hole pair associated with the tt 
bands at the K-point. This is in reasonable agreement with the values obtained from other 
experimental probes. As a corollary, this model can also explain the difference recently 
observed in the peak positions of the D and D* bands in Stokes and anti-Stokes spectra. 
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FIGURES 

FIG. 1. D-mode frequency dependence on excitation energy for different forms of car- 
bon-disordered carboni (filled squares) , glassy carbon^ (filled triangles) and carbon nanotubesi 
(open circles). 

FIG. 2. Raman spectra of graphite at E^, = 4.82 and 2.54 eV. The inset shows the intensity 
ratio of D-band to that of G-band as a function of as reported by Wang et aZ.i. The peaks 
marked by asteriks are the internal vibrational modes of sulphate ion for use as internal intensity 
calibrant. 

FIG. 3. (a) Schematic band structure of graphite showing the four vr bands at the K-point of 
the Brillouin zone. The wavevector k is perpendicular to the c-axis. (b) Schematic of the single 
graphene sheet at the K-point. Dashed lines denote the Fermi energy. 

FIG. 4. Phonon dispersion curves in graphite calculated along certain high symmetry axes as 
given by Mathew et aZ.0. 
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